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It has been suggested that superhydrophobic surfaces, due to the presence of a no-shear zone, can 
greatly enhance transport of surface charges, leading to a considerable increase in the streaming 
potential. This could find potential use in micro-energy harvesting devices. In this paper, we 
show using analytical and numerical methods, that when a streaming potential is generated in 
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^vq ' such superhydrophobic geometries, the reverse electro-osmotic flow and hence current generated 

by this, is significant. A decrease in streaming potential compared to what was earlier predicted 
is expected. We also show that, due to the electro-osmotic streaming-current, a saturation in 
both the power extracted and efficiency of energy conversion is achieved in such systems for 
large values of the free surface charge densities. Nevertheless, under realistic conditions, such 
microstructured devices with superhydrophobic surfaces have the potential to even reach energy 
conversion efficiencies only achieved in nanostructured devices so far. 

V 

i 

Energy conversion and flow detection devices, fabricated in the micro- and nano-scale, have a 
potential application in novel highly integrated or portable systems. One way to achieve conver- 
sion from hydraulic pressure to electric energy is via the streaming current. The phenomenon 
*£h was first reported by Quincke 1 in 1859, who noted that when pure water flows trough a porous 

material an electric current is generated. Helmholtz-, expanding on Quincke's reasoning, gave 
the first theoretical explanation for this, noting that virtually any surface brought in contact 
with water will acquire a net charge, attached to the surface, that is screened in a diffuse charge 
cloud within the water in the vicinity of the surface. Thus when fluid flows through a capillary, 
£ — , the mobile part of this electrical double layer is dragged along, leading to the streaming current. 

However, since the electric double layer naturally arises at solid walls where the flow velocity 
, is small, conversion rates are usually low 3 - - — , reaching a few percent in channels of nanome- 

ter dimension-^ and much less in microchannels. It was therefore suggested that an effective 
wall-slip can increase the efficiency of such a device^— . 

Some superhydrophobic surfaces trap pockets of air within grooves. This so-called Cassie- 
Baxter state has attracted much attention due to the reduced drag experienced by a fluid flowing 
over such a surfac e - 13 ' 14 . In particular shear and pressure driven flow over a striped geometry 
has been thoroughly studied^ - — . Since the air-water interface can be considered as posing no 
resistance to shear, in principle such a surface seems ideal for an increase in streaming current, 
as long as the free surface carries a net charge. This latter condition seems nontrivial at first, 
since any ions that are preferentially attracted to the surface are expected to attract an equal 
amount of opposite charge in a diffuse cloud just as on a solid surface. However, a net charge is 
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attracted to the surface when an electric field is applied normal to the fluid-air interface**, for 
example by embedding electrodes within the grooves 21 . 

It has been shown that the super hydrophobic nature of the surface, together with a charged 
air-water interface, can be exploited in order to generate a significant increase in streaming 
potential^. However, in this study the effect of electro-osmotic (EO) flow due to the built-up 
streaming potential was neglected, which generally decreases the efficiency since it goes in the 
opposite direction as the pressure driven flow—. Moreover, on superhydrophobic surfaces the EO 
flow is enhanced by essentially the same mechanism and magnitude as the streaming current, 
namely due to charges on the no-shear surface 2 ^ - — . In this paper therefore, we show that 
on superhydrophobic surfaces, the electrokinetic flow resulting from the streaming potential, 
together with its streaming current, plays an important role in limiting the efficiency of energy 
harvesting with increasing free surface charge and effective slip length. So far, the only paper 
where both effects, in particular the EO streaming current, have been properly taken into 
account for energy conversion on patterned surfaces seems to be Ng and Chu's work 2 ^. Their 
focus, however, is on the EO flow and the streaming current due to pressure driven flow, only 
briefly exploring energy conversion on a patterned surface with alternating no-slip and finite- 
slip stripes. Here, our attention is on energy conversion in the optimal case of alternating slip 
and no-slip surfaces, a good approximation to a liquid in Cassie-Baxter state above air-filled 
cavities. Further, in our analysis we numerically go beyond the Debye-Hiickel approximation 
and explore its impact on the results for energy conversion. 

We organize this paper into three sections. In section [H] we present the general formulations 
and derive expressions for the pressure driven flow, the EO flow and energy conversion and 
harvesting. In section IIIIl we present the results of our study in two parts. We first discuss 
our results through analytical approximations of the expressions derived in section HH and then 
graphically present the exact behavior of the expressions presented in the section [Hi and validate 
our analytical results. In section TV\ we present our conclusions based on this study. 

II. GENERAL FORMULATIONS 

We consider flow over a striped surface with alternating no-shear and no-slip regions as shown 
in figure[TJ The liquid-air interface is assumed to be ideally flat, such that the liquid occupies the 
domain y > 0. Zhac 5 showed that the overall efficiency for generation of streaming potential 
is greater for flows in the longitudinal direction when compared to the transverse direction. 
Therefore, we here concentrate on flow in the longitudinal direction. The electric potential on 
the no-slip part of the superhydrophobic surface as well as on the no-shear part is assumed to be 
constant, denoted by ( ns and ( s , respectively. The striped surface is assumed to have a period of 
w and the height of the channel is taken to be 2h, with the channel being symmetric about the 
center line. Since the device is assumed to be in the urn range, and the typical Debye screening 
length is of the order of nm, we make the assumption that (h S> \d) where Ad is the Debye 
layer thickness and n = 1/Xd- This condition will be relaxed later in numerical calculations. 
The liquid viscosity is denoted by \x and the permittivity of the liquid is denoted by e. The 
free surface fraction, i.e. the ratio of the width of the no-shear region to the total width of the 
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FIG. 1. Schematic diagram of the superhydrophobic surface. Within the grooves air is entrapped and we assume 
the air-liquid interface to remain ideally flat. The length of the channel is assumed to be I and consists of two 
parallel surfaces at distance 2ft. The surfaces are periodically structured with groves of width w8 and periodicity 
w. 

period of the interface, is denoted by 5. A schematic of the setup is shown in figure [TJ 

A. Bi-Layer Potential Distribution 

For our (semi-) analytic calculation, the potential distribution in the bi-layer is derived based 
on the Debye-Hiickel theory^. The electric potential in the system obeys the Poisson-Boltzmann 
distribution. Under the assumption that the electric potential energy of the system is much less 
than the average thermal energy, ksT/e ~ 26 mV at room temperature, the governing equation 
for the potential distribution is 

V 2 V(x,y) = K 2 V(x,y). (1) 

This equation is solved under the boundary conditions that, V(x, 0) = Cns for < x < ^ and 
V(x, 0) = Cs for < x < ^, using separation of variables and a Fourier series expansion in the 
periodic variable x. Denoting the average potential on the boundary as ( = 5( s + (1 — 5)( ns , 
the final form of the solution is given by 

oo 

V(x, y) = Ce~ Ky + Pn cos (u n x) exp y) , (2) 

n=l 

for y < h, where we have used uo n = 2 ^ and (3 n = 2 ^ a ~^|_ hm " 7r ^ i Here and in the following, we 
will assume nh S> 1, such that exponential functions ~ exp(— ny) have died away sufficiently in 
the center of the channel, allowing us to treat the walls independently for all solutions of this 
type. 

Later, we also use a simplified form of equation ([2]), obtained by solving equation ([1]) inde- 
pendently in the no-shear region and the no-slip region. The resulting solution is defined in the 
form of a split function as 



V(x, y) = Cse~ Ky , < x < 5w/2, 
V{x, y) = ( ns e- Ky , 5w/2 <x< w/2. 



(3) 
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This approximates equation ([2]) in a region extending O(Xd) from the wall. We will comment 
on the applicability of this approximation in section Hi CI where it is used. 

The charge-density corresponding to the potential distribution is frequently used and is given 
by Gauss's law 



Pe 



-eV 2 F = -ek 2 V. (4) 



B. Streaming Velocity and Current 



The streaming current is generated by pressure-driven flow. Since flow velocities are small, 
inertial terms can be neglected in the momentum equation and the governing equation for the 
streaming velocity thus is the Stokes equation 

AP 

V 2 Up (x,y) = — , (5) 

where AP is the (negative) pressure drop applied across the channel of length /. This equation 
is solved for the velocity profile by splitting the velocity into two components, u p = uL + u 2 . 
Here, uj, is the velocity profile for flow between parallel no-slip plates. The correction, u 2 , to 
this velocity is then obtained by demanding that u p obeys the no shear boundary condition on 
the slipping portion. The flat plate velocity profile is given by 



u p\v) = 



2/zZ 



2 y__(y_ 

h \h 



(6) 



The correction to this velocity profile satisfies V 2 u 2 = 0. Since the geometry is periodic, the 
velocity profile can be expanded in a Fourier series as 

oo 

Up(x, y) = a + a n cos(^ ra a:) e~ w " y . (7) 

The net velocity profile, u p , satisfies the boundary conditions that at y = 0, d y u p = in the 
no shear region and u = in the no slip region. The velocity correction, ui, thus satisfies the 
conditions at y = that d y Up = —d y up; < x < ^f- and n 2 = 0; < x < f . Using the 
expansion ([7]) these equations translate to 



^ APh Sw 

> u n a n cos uj n x = — , < x < — , 

i fd 2 

n=l 



oo 



CLQ + ^2 a n COS UJ n X = 0, 



n=l 



5w W 

< X < — . 

2 2 



(8) 



Systems of equations of this type were solved by Sneddon^!, such that^ 

APwh 

APwh f Sn x 



2w(J,l 



f x 

/ tan— [P n (cos x) + P n _i(cos x)\ dx, 
Jo 2 



(9) 



where P n (x) are the n th order Legendre polynomials and /3ii = ^ log sec (Jy) is the non- 
dimensionalised effective slip length for pressure driven or Couette flow over such surface o 15116 . 
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The electric streaming current through a cross section spanning half a period generated by 
the pressure driven flow is given by 

fh 

p e Updxdy. (10) 



Jp — 2 



o Jo 



Using the expressions derived in equations 
give 

eK 2 APhw( ~ 
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and ([9]), the above integral is evaluated to 

a n f3 n w 



2 00 

-Y- 

2 n=0 Un + V* 2 +Un 



(11) 



Since the typical order of Ad = 1/nis much less than the width of the plate, we can approximate 
the sum in the above expression by 



~0~ 



a n (3 n w 



n=Q UJ n + ^K 2 + £j2 



ek(C s - ( ns )w ^ sin(n7T(5) 
" 2 > Q "' 



AP 



n=l 



717T 



(12) 



By integrating equation ([8|) from ^ to ^, we can show that Ylri=i a " S ' D lT^ = (1 — ^) a 0- 
Substituting this in equation (TTZj) we get a simplified expression for the streaming current 



J p ~ w/i 



AP. 



2/j,l pi 

Correspondingly, the volumetric flow rate for the pressure driven flow is obtained as 



Q P = 2 



h 3 wAP (2 
u p dx dy = - + a ) ■ 

o Jo V V 3 



(13) 



(14) 



C. Electro-Osmotic Flow and Current 



An electro-osmotic flow, which thereby drives a current, results from the electric field which 
is generated due to the pressure driven flow. The governing equation for this EO flow is given 
by 



V 2 u e = p e E z . 



(15) 



We again write the velocity profile, u e = u\ + u 2 , as a superposition of a velocity on a no-slip 
surface with periodic potential, u\, and a correction term which arises because of a no-shear 
zone, tig. Assuming that d 2 u], -C d 2 ul, correct apart from a small region of width of order Ad 
around x = Sw/2, we solve for u\ using equation ([3]) for the electric potential and get 



ul(x,y) = £ -^(l-e-y) 



< x < 



5w 
2~' 



5w w 
-2 <X< r 



(16) 



Obviously, this is not a full solution of ([T5]) . since it has a jump at x = Sw/2, just as ([3]) is not 
a full solution of (HJ. However, as we will see, the EO velocity field is generally dominated by 
u 2 when charge is present at the free surface, so the error introduced due to this approximation 
is exceedingly small. Moreover, and more importantly, for thin Debye layers, kw 3> 1, the 
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approximation is excellent in the vicinity of the wall, such that we can obtain the correct 
boundary condition for u 2 from it. Effectively, this amounts to replacing the diffuse charge 
cloud with a corresponding surface charge density at the interface, resulting in a jump in shear 
rate across the interface^. 
The general form of u 2 can again be expanded as a Fourier series 



(x,y) = A G + ^A n cos{oj n x)e-^ y . (17) 



n=l 



Since the full EO velocity field, u e , satisfies d y u e = at y = in the no shear region and u e = 
in the no slip region, the boundary conditions for u 2 at y = are d y u 2 = —d y u],; < x < ^ 
and u 2 = 0; 4*r < x < In particular, the Fourier coefficients of u 2 must satisfy the conditions 



} LU n A n cos u n x = — , < x < — , 

i fil 2 

n=l 

OO r- 

x - OW W 

A + 2_^A n cosu n x = 0, — < x < — , 

n=l 

completely analogous to equation Q. Solving this double series, we get 

£( s KwA(j) 



(18) 



. £( s KwA(j) f I 
^ = n i / taI1 n "n(cOS X) + -P n -1 (COS X) dX. 

2-Kpd J 2 



(19) 



We now see that the homogeneous part of the flow, given by Aq, is larger than u], by a factor of 
kw[3\\/2. Thus as long as ( S /Cns *C /3\\w/\d, the approximation going into equation (fT6|) plays 
no role in the bulk velocity field. 
The electro-osmotic current, J eom , is thus given by 

Jeom = 2 / p e u e dx dy. (20) 
J JO 

This integral is evaluated in Appendix El and the expression for the EO current is given by 
equation ()A2|) and is shown below 



2\ 



egA^w + eC^A0(l - 5)w + C^£CgA0 2 ^ ^ 



4«£t£ 4k/X£ AfiL 7T \ 2 J 

■m 

Under the assumption that (k 3> w n ), a simplification to the above expression is carried out, 
similar to the simplification of the streaming current, and is given in Appendix^] The simplified 
form of the EO current is given by 

'e 2 ( 2 K 2 w 2 ' 



Jr., 



Acj). (22) 



Correspondingly, the flow rate due to the EO velocity field is given by 

Q e = 2 [ 1 2 u e dxdy = -wh ( A + A _ J_ ] j . ( 2 :>>) 



JO 
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(25) 



D. Energy Efficiency and Streaming Potential 

In our system, apart from the streaming current and the electro-osmotic current, a bulk con- 
duction current through the system is present which is given on the basis of the bulk conductivity 
of the liquid, a, by 

J cond — ^ j Kr*) 

valid for homogeneous conductivity, i.e. neglecting variations in carrier density within the double 
layer, a condition we will relax in our numerical evaluation. Energy is extracted from the system 
by connecting an external load of resistance R ex t in series with the channel. The current running 
through this load is thus J ex t = A(j)/R ext and the power extracted is P ex t = A(j) 2 /R ext . 

Due to the linearity of the system, the expressions for the currents and flow rates in the 
system, derived in sections III Bl and III CI can be succinctly represented in matrix form as 

( Jnet\ = (L 4> L p \ ( A A 
\Qnet) S p ) \AP) ' 

In the above representation, = L eom + L con d, is a sum of the contributions due to the EO 
current and the bulk conduction current. The coefficients given in equation (|25p have been 
calculated in equations (fTT|) . (|2~Tj) . ([M]L flT3J) and ([23]) . Based on general considerations in non- 
equilibrium thermodynamics it was show n 28 i 29 that the matrix in the above representation is 
symmetric, i.e. L p = S^, a relation referred to as Onsager reciprocity. We will later see, that 
this is in fact valid in our approximation and was also obeyed in our numerical simulations. 
Note that this Onsager relation has been elegantly shown to be valid on striped surfaces by Ng 
and Chi>2£ and as such can be used in checking the validity of numerical calculations. 

Kirchoff 's law dictates that the currents generated in the system be balanced. Since the sign 
of the EO current is negative, we have 

L p AP+^ + (L cond + L eom )A0 = O. (26) 

The streaming potential is given by A<j) when R ex t oo. Conversely, when we want to draw 
maximum power from the circuit, R ex t must satisfy R~ xt = L eom + L con ^. Thus, the maximum 
power which can be extracted from the circuit is given by 

Pmax = — ~ = (Atp) 2 (L con d + L eom ). (27) 
i^ext 

The input power into the system is derived from the pressure-driven flow as Pi n = AP(Q P +Q e ). 
The efficiency of energy conversion, r/, is the ratio between harvested and input power, rj = 
Pout /Pin, and is dependent on the external resistance used for extraction. Based on the Onsager 
reciprocity, the efficiency can be elegantly written in terms of the figure of merits, Z, 



Z = (28) 



and the maximum efficiency then becomes 



Z 

^ ~ (1 + ^1— Z)2- ( 29 ) 
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III. RESULTS AND DISCUSSION 

Based on the earlier formulations, the effect of electro-osmotic flow on the energy conversion 
efficiency, output power and the streaming potential developed were calculated. The results are 
presented in two parts: first we look at the Onsager relation and approximate analytical results; 
then we look at a more exact graphical representation of the solution. 

A. Approximate analytical solutions 

The Onsager relation predicts that the amount of current generated per unit pressure drop, 
L p , should be the same as the flow generated per unit voltage drop for the system, Sj,- The 
coefficients corresponding to the current generated by the pressure field and the flow generated 
by the potential are given by equations (fT3|) and ([23]) respectively. Thus, we see that under the 
thin bi-layer assumption (nh 3> 1) 



Lp = S<j> = wh 



eCsKW s( 



(30) 



2fj,l 11 ill 

in accordance with the Onsager relation. In the same limit, the diagonal coefficients in equation 
([25]) are given by 

L eom =( £ -^^)p h (31) 



2{il 

L C ond = — I — • (32) 

For reasonably large values of the free surface fraction one has kw/3\\ 3> 1, so the first term in 
equation (|3U|) will dominate this expression and the last term may safely be neglected, which 
we will do in the following. Moreover, instead of using the potentials £ s and ( ns on the free and 
no-slip surfaces, we will mostly refer to the corresponding surface charge densities, a s and a ns , 
which in the Debye-Hiickel approximation are related by 

(J s = sQgK,, <T ns = E^ ns K. (33) 

Substituting these quantities in equations ([27]) and fl2SD, we obtain simplified analytical expres- 
sions for the maximum power and figure of merit, Z, 

AP 2 w 2 h 2 p l{ 1 

Pout ~ 8fjH 7] | o/^M ' ' 1 1 ! 



(35) 



From the expression for the figure of merit, Z, the efficiency of power extraction can be easily 
obtained using equation (|29p . In our approximation, since transport on the free surface is 
completely dominating the system, both the maximum power extracted and the efficiency of 
extraction are independent of the charge on the no-slip surface and the Debye layer thickness. 
For large free-surface charges, a saturation behavior is observed, c.f. equations (fM|) and (|35|) . 
This occurs because of the linear dependence of the streaming current on the charge distribution 
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or the free surface charge, while a quadratic dependence exists in the relationship of the EO 
current on the free surface charge, since the EO flow which drives the current is itself dependent 
on the surface charge. The efficiency at saturation depends only on the ratio of (effective) 
length scales, h/(w/3\\), and its dependence on the free surface fraction 5 is graphically seen 
in section fill CI In particular, it only depends on a geometric ratio of (effective) length scales, 
h/ (w/3h), and the dimensionless group /iah/ (of w^A, which is the ratio between charge transport 
by bulk conduction and electro-osmotic streaming and determines how close the efficiency (|29p 
and power can come to their saturation value. 

It is a simple matter to calculate the corresponding efficiency for a parallel plate arrangement, 
where instead of the no slip condition a Navier slip condition, bd y u = u, with slip length b is 
used at the walls. In the limit k > 1 and 6k > 1, the results in this case are completely 
analogous to equations flM]) and ([35]) when we make the identification b = f3\\w/2. Note that 
this is in complete agreement with and directly related to the finding, that for a striped surface 
b = (3\\w/2 acts as an effective slip length. 

Here we again like to stress the fact that these results, equations (|34|) and (|35|) . are valid in the 
limit of small Debye lengths, Ad <C w, for not too large ratios of ( S /Cns "C /3\\w/Xd and only for 
h> w. Although the result was obtained in the Debye-Hiickel approximation, ( < ksT/e, this 
assumption is likely less stringent, as long as we interpret eCsk as the surface charge density in 
all occurrences, since this quantity is more directly related to the driving force at the interface. 
In section IIIICl we will compare these results with the more exact ones based on expressions 
(fTTj) . (HU), ([2"T]) and ([2"3"]l . Additionally we also compare with numerical calculations, using 
the nonlinear Poisson-Boltzmann equation instead of relying on the linearised Debye-Hiickel 
approximation . 

B. Numerical Methods 

In order to probe the range of validity of our (semi-)analytical results, we use a commercial 
finite element package, Comsol Multiphysics^, to go beyond the linear Debye-Hiickel ap- 
proximation, eq. (H|). In particular, the electric potential is obtained by solving the nonlinear 
Poisson-Boltzmann equation appropriate for a symmetric electrolyte with ions of equal mobility, 



where Vq = ksT/e ~ 26 mV. On the no-slip and no-shear boundaries the charge density is 
prescribed, corresponding to a Neumann boundary condition 



We use here the charge density instead of the zeta potential, since it more directly reflects the 
driving force for the EO flow on the free surface. We remark that the appropriate conversion 
between charge density and zeta potential at large potentials can, instead of equation ([33]) in 
the linear regime, be approximately achieved via the Gouy-Chapman solution available for a 
wall with fixed charge density^. 




(36) 



ed y V\ y= Q = CTj, (Ti 



a s , < x < Sw/2, 
a ns , 5w/2 < x < w/2. 



(37) 
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The governing equations for pressure driven and EO flow remain the same, equations ([5]) 
and (|15p . with p e = — sk 2 Vo sinh(VyVo) as charge density, however. As before, the streaming 
current is obtained by integrating the product of velocity field and charge density, equations 
()10p and (|20p . Similarly, ion concentrations vary according to the Boltzmann equation, Cj = 
Cbulk exp(— ZiV/VO), where c^uik is the bulk ion concentration and z% = ±1 is its valency. Hence, 
the local conductivity is given by a(x,y) = abulk cosh(y/V0), where abulk = HECbuik with 
ion mobility fiE and Faraday's constant F. Practically, unless stated otherwise, we prescribe 
certain values of Ad and abulk and from these the corresponding values of bulk concentration, 
Cbulk = ^^b/(2i ? A^ ) ), and mobility, \ie = Gbulk^b I ( e ^b) > are inferred. The total conduction 
current thus becomes 



replacing equation (24), which is used apart from cases where we explicitly want to highlight 
the influence of the Debye layer conductivity. Thus all terms of the transfer matrix (|25p are 
numerically known and the efficiency analysis can proceed as before. 

A typical computational grid consists of N x x N y = 100 x 300 elements. Particular care was 
taken to refine the grid close to boundaries at y = in order to capture the steep gradient 
in electric potential and EO velocity and obtain grid independent results^. A fifth order 
discretization was used for all the PDE's. 

C. Graphical Analysis 

In this section, we plot the behavior of the analytical expressions for maximal power and 
efficiency derived in section IIII Al in some illustrative examples. These results are compared 
both with the more exact semi-analytic expressions presented in sections III Bl III CI and III Dl 
as well as with numerical data to asses the influence of a finite Debye-layer thickness. We 
assume the fluid properties to be those of water, i.e. [i = 1 mPas, e = 80 • 8.85 10~ 12 As/(Vm). 
Initially, we will consider a low bulk conductivity of a = 0.1 (J.S/cm of highly purified water 
and thus also a correspondingly large Debye length of Ad = 1 (J.m for the comparison. Note 
that for our numerical implementation this corresponds to the hypothetical mobility of [Xe — 
5.43 • 10 m 2 /(Vs) and bulk concentration Cb u ik — 10 -7 mol/L. 

Figure [2] shows the variation of the output power with amount of charge on the free surface. 
With increasing surface charge densities on the free surface the power extracted increases until 
the electro-osmotic effect begins to dominate. Subsequently, a saturation in the power extracted 
is achieved. The power output at saturation and the surface charge at which the saturation is 
reached varies depending on the amount of free surface present, parametrised by the free surface 
fraction 5. The dashed line in this figure further illustrates the role played by the electro-osmotic 
flow and its streaming current: when EO flow is neglected, the power extracted from the circuit 
grows rapidly and without bound with increasing amount of charge on the free surface. Further, 
our semi-analytic and numerical data suggest that the analytical solution accurately predicts the 
behavior of the system given a constant in bulk conductivity of the system. In order to highlight 
the role played by the increase in conductivity within the Debye layer, we have evaluated the 




(38) 
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FIG. 2. Output power density for varying free surface charge density, a s , for different values of the free surface 
fraction, S. The solid line is calculated based on our analytical expression (I34p . and the crosses (x) are calculated 
based on the more exact semi-analytic expressions (|14l) . (12ip and (|23|) . Pluses (+) and circles (o) are 

numerical results using the nonlinear Poisson-Boltzmann equation, differing in the evaluation of the conduction 
current by, respectively, a constant and variable charge density according to equations (|24p and (|38[) . The dotted 
line shows the behavior when the electro-osmotic flow is not accounted for (at 8 = 0.7). The figure was plotted 
assuming w = h = 100 u.m, a — 0.1 pS/cm, Xd = 1 M- m > /x = 1 mPas and a s /a ns = 2. 



conductive contribution to in the numerical data both assuming a constant charge density, 
equation (|24p . as assumed for the (semi-) analytic evaluation, as well as with a variable charge 
density according to equation (f38jh It is evident that here there is little difference between the 
two cases. 

The role played by the EO flow and its streaming current is further illustrated in figure O 
where the relative drop in streaming potential compared to the case without EO flow is shown. 
We find that the drop in streaming potential when we account for the EO flow is quite significant, 
especially for higher values of the free surface charge. Further, this drop increases for increasing 
ratios of the free surface fraction, 5. This increasing contribution of the EO current results in 
decreasing the streaming potential observed for increasing surface charge. Note that in figure 
[3] compared to figure [2] we consider a smaller range of charge densities at the free surface and 
that here the difference between analytic and semi-analytic results is negligible. 

Figure H] shows the variation of efficiency for two different values of free surface fraction 5 
and for different values of bulk conductivity a. We see that the maximum efficiency which can 
be achieved for these devices is independent of the bulk conductivity of the fluid, as evident 
from equation (|35p and (|29p : in particular, it is dependent only on the geometry of the channel, 
such as the width, height, length and the free surface fraction. However, we see that for lower 
values of a, the saturation efficiency is achieved for lower values of the charge density on the 
free surface. We note that at fixed surface charge density, the efficiency-increase with 5 is not 
linear and grows more rapidly the closer 5 becomes to 1; a similar trend is followed by the power 
output, as shown in figure Again the analytic expression agrees well with the numerical data. 

We would like to stress here that, when attracting charges to the free surface by a normal 
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FIG. 3. Relative drop of streaming potential for increasing free surface charge density, a s , for different values of 
the free surface fraction, S. The solid line is calculated based on our analytical expression (|34[) . and the crosses 
(x) are calculated based on the more exact expressions |TT}, I] 14ft. (|21l) and (|23[) . The figure was plotted assuming 
w = h = 100 (xm, a = 0.1 u.S/cm, Ad = 1 um, = 1 mPas and a a /a ns = 2. 
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FIG. 4. Variation of efficiency with free surface charge density, a s , both for varying values of bulk conductivity 
a and free surface fractions S = 0.4 and S — 0.9. Lines are calculated based on the analytical expressions (|29l 
[35}. Circles (o) are numerical results obtained using the nonlinear Poisson-Boltzmann equation. The dashed 
lines shows the maximum theoretical efficiency predicted for such systems. Here w = h — 100 urn, A_d = 1 urn, 
fj, = 1 mPas and a s /cr ns = 2. 



electric field as described in the introduction, a charge density of 10 -4 C/m 2 converts to an 
applied field of E ext = a s /eo m 10 7 V/m, about three times the breakdown field strength for 
dry air. To reach such high surface charge densities by this method, it would thus be necessary 
to use a medium with higher breakdown voltage or larger permittivity e. (For a discussion of 
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Charge Density on the Free Surface (C/m 2 ) x 1 ry 4 D V h 

(a) (b) 

FIG. 5. (a) Variation of maximum efficiency with free surface charge density, cr a , for values of \o/h between 
0.01 and 0.1, fixing the geometry length scale h and varying A_d. Lines are calculated based on the analytical 
expressions (|29l I35[) . Circles (o) are numerical results obtained using the nonlinear Poisson-Boltzmann equation. 
Here w — h = 100 urn, 8 = 0.9, fis = 5.43- 10~ 7 m 2 /(Vs), = 1 mPas and a B /a na = 2. (b) Variation of maximum 
efficiency as a function of Xofh at fixed Au = 1 um, a — 0.1 uS/cm and variable length-scale h for several values 
of surface charge density a a . Circles (•) are numerical results obtained using the nonlinear Poisson-Boltzmann 
equation and the dotted line is the analytical expression I35|) . Here w — h, 8 — 0.9, fi — 1 mPas and 
u s /a ns = 2. 



other electric breakdown mechanisms on super hydrophobic surfaces we refer to Oh et al.— .) 
Nevertheless, for large enough free surface fractions, efficiencies of the order of a percent seem 
possible for a microstructured system by this method, something that otherwise is only achieved 
in nanostructured devices^. 

The graphs presented above represent a case where the periodicity of the surface pattern and 
the height of the channel were assumed to be 100 um, and the fluid properties are based on 
those of highly purified water. A low conductivity of the system was taken in all cases. In figures 
[2] and El the Debye layer thickness was assumed to be a high value of lixm, corresponding to the 
low concentrations typically present in low conductivity aqueous systems. To a large degree, 
the behavior of the curves is independent of the Debye Layer thickness and the charge on the 
no-slip surface, since our analytical expressions conform with the predictions based on our more 
exact formulations. This slight dependence on the Debye layer thickness disappears for smaller 
values of and our analytic results, (fMj) and (f35|) . become excellent approximations. We find 
that the efficiency is determined mainly by the ratio of f3»w/h. Our results thus indicate that 
for microfluidic devices with increasingly large portions of free-surface, we would be able to 
achieve much higher energy conversions due to the steep increase in effective slip length— . 

To gain further insight into the dependence of the efficiency of energy conversion for increasing 
thickness of the Debye-layer additional numerical calculations were performed. For this we fix 
the ion mobility, /ig, and prescribe variable values for Ap, in each case choosing the bulk 
concentration and corresponding conductivity in agreement with these values. In this way we 
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study a fixed electrolyte at varying bulk concentration (and thus varying bulk conductivity). 
Since we are mainly interested in the impact of Ad on our results and less in a particular 
electrolyte system we fix the mobility at the hypothetical value of /jle = 5.43 • 10 _7 m 2 /(Vs), as 
used previously. The results of these calculations are found in figure [5] (a) for values of Xd/w 
varying from 0.0075 to 0.05 and are compared with the analytical expression, ([34"]) and ([35]) . At 
fixed surface charge density the efficiency of energy conversion increases for increasing values 
of the Debye layer thickness. This just reflects the fact that an increasing Xd corresponds 
to decreasing ion concentrations and thus conductivity. As before the values saturate at the 
the value predicted by our theoretical expressions for large values of the free surface charge. 
In addition to that, even below saturation, the analytic expression, equations ([35]) and ([29]) . 
approximates the numerical results well. However, it is apparent, that the deviation will become 
larger with a further increase in Debye length. Nevertheless, this indicates that the validity of 
our analysis also extends to aqueous systems in smaller structures than presently considered. 

To demonstrate this more directly, we also vary h at fixed Xd and bulk conductivity. This is 
shown in figure [S] (b) for the case of highly purified water, i.e. Xd = 1 [im and a = 0.1 U-S/cm, 
and for several values of surface charge density a s between 2.5- 10~ 5 and 5- 10 _4 C/m 2 . For small 
values of the free surface charge the graphs shows the characteristic"^ increase of the efficiency 
with Xnjh to a maximum and the decrease of efficiency upon further reduction of the length 
scale due to increased average conductivity in the channel. For increasing a s the maximum is 
shifted towards smaller values of Xo/h, until it vanishes and a monotonic decrease is observed 
for a s > 7.5T0 _5 C/m 2 in this case. Together with figure[5](a) this gives a fairly good indication 
for the range of applicability of the analytic expression, ([34"]) and ([35]) . Nevertheless, since our 
main interest lies with larger structures we will not pursue this matter further here. 

IV. CONCLUSION 

In this paper we have analyzed the energy conversion from hydraulic pressure to electric 
energy in a device with large free surface fractions, as observed for a liquid in Cassie-Baxter 
state on superhydrophobic surfaces. Due to the substantially reduced drag, pressure driven 
flow is greatly enhanced over such surfaces. Accordingly, when charges are present on these 
free-slip surfaces one obtains a greatly enhanced streaming current, proportional to the surface 
charge density. However, in this case also the EO flow and its corresponding streaming current 
is enhanced by the same mechanism, depending linearly and quadratic on the charge density 
on the free surface, respectively. Due to this quadratic dependence, the EO streaming current 
strongly influences and limits the efficiency of the device, leading to a saturation of efficiency 
for large charge densities on the free surface. The cross-over when EO flow becomes important 
is dictated by the dimensionless group fiah/(agW/3^), the ratio between bulk and EO streaming 
current. As we have shown, this is equivalent to a device with unstructured surfaces obeying a 
Navier-slip condition with a slip length of b = f3\\w/2. 

This observed saturation with increasing charge density on the free-surface implies that it 
would be of greater importance to design devices with a higher free surface fraction, which 
would lead to a considerable increase in the efficiency of energy conversion. These microstruc- 
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tured devices have the potential to be far more efficient than flat plates with no slip present 
everywhere and can even reach efficiencies only achieved in nanostructured devices so far. We 
await further experimental evidence to support the use of such systems both for energy con- 
version in microfluidic chips and as flow-rate detection devices integrated in conduits or even 
made for implants on surfaces such as blood vessels. 
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Appendix A: Simplifying the Electro-Osmotic Current 



The EO streaming current, equation (|20p . can be simplified along the same lines as the 
streaming current due to pressure driven flow in equations (|llj) - (|13|) . Using the definition of 
the potential distribution provided in equation Q and the expressions for the electro-osmotic 
velocity given by equations (fT6j) and (fT7|) . we simplify equation ([20]) to get 
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Performing the integrations in equation (|A1|) . leads to 
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Assuming that kw 3> 1, we can approximate equation (|A2j) as 
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By integrating equation (fTgj) from to 4^, we can show that Y^=i 
Therefore, equation (|A3|) can be written as 
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For reasonably large values of S where log sec 4^ ~ O(l), we see that because wk > 1 we have 
log sec 4^ S> d + (^j) 2 (l — d). Therefore, J eom can be approximately written as 
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